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COMPUTATION  OF  B£AM  PATTERNS  AND  DIRECTIVITY 
INDICES  FOR  THREE-DIMENSIONAL  ARRAYS  WITH 
ARBITRARY  ELEMENT  SPAC1NGS 

INTRODUCTION 

A  special  formulation  of  the  computation  of  beam  pattern  functions  of 
arbitrary  arrays  is  presented  herein,  along  with  3ome  computational  aspects. 
A  package  of  computer  programs  has  been  written  in  FORTRAN  V  and  is  de¬ 
scribed.  Complete  input  parameter  definitions,  output  information,  usage,  and 
pregram  descriptions  are  given  in  detail  in  the  text  and  in  the  tables.  Selected 
numerical  examples  exhibit  the  different  options  and  show  the  typical  results 
obtained  with  the  beam  patterns  being  described  graphically  and  the  directivity 
index  against  noise  and  reverberation  being  given  numerically  in  decibels.  The 
computations  associated  with  the  program  are  carried  out  on  the  Univac  1108 
computer  in  single-precision  arithmetic  with  the  EXEC  8  system.  A  listing  of 
the  computer  programs  is  included  in  the  appendix. 

The  package  contains  some  important  features  which  are  worthy  of  note. 
First,  the  design  allows  treatment  of  arrays  whose  elements  are  distributed  in 
three-dimensional  space  in  a  completely  arbitrary  fashion,  with  either  uniform 
or  nonuniform  spacing  between  elements.  Determination  of  beam  pattern  func¬ 
tions  is  computed  in  a  maimer  that  allows  the  beam  to  be  defined  in  any  plane 
surface  which  passes  through  the  array  center.  By  a  selective  sectioning  in 
angular  spherical  coordinates,  this  technique  can  produce  representations  cf 
three-dimensional  beam  shapes. 

The  program  provides  for  array  response  to  either  plane  or  spherical  wave 
excitation.  In  addition,  in  order  to  treat  completely  arbitrary  arrays,  an  indi¬ 
vidual  array  element  response  is  incorporated  which  has  directional  receiving 
properties.  Further,  the  program  is  capable  of  incorporating  any  amplitude 
shading  'cciJdque  desired.  Arbitrary  phase  delays  in  element  response  may  also 
be  added,  apart  from  the  usual  artificially  introduced  array  steering  delays. 
This  auditions-.!  feature  is  useful,  for  example,  in  determining  array  perform¬ 
ance  when  unwanted  or  uncompensated  phase  delays  are  present  in  the  system. 

Directivity  index  and  reverberation  index  are  computed  from  beam  pattern 
function 3  by  the  application  of  numerical  quadrature  techniques.  An  option  is 
prodded  for  integration  in  subintervals  with  different  mesh  sizes  in  order  to 
minimize  computer  running  time. 
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MATHEMATICAL  FORMULATION 


BEAM  PATTERN  FUNCTION 

The  beam  pattern  function  is  considered  to  be  dependent  on  the  two 

spherical  coordinate  deviation  angles:  $,  defining  the  azimuthal  deviation  in  the 
horizontal  plane,,  and  4>,  the  vertical  deviation  from  a  horizontal  reference.  The 
reference  direction  in  azimuth  may  generally  be  assigned  to  coincide  with  a  geo¬ 
metric  axis  of  the  array  if,  as  is  usually  the  case,  there  is  inherent  symmetry 
in  the  physical  array  structure.  The  origin  of  coordinates  is  assumed  to  be  the 
geometric  center  of  the  set  of  elements  comprising  the  array,  as  shown  in  fig¬ 
ure  1.  The  array  element  is  designated  by  its  Cartesian  position  (x^,  y^,  z^) 
k  =  1,  2,  . . . ,  N. 

ARRAY  RESPONSE  TO  PLANE  WAVE 
EXCITATION 

A  plane  wave  that  is  impinging  mi  the  array  is  completely  defined,  insofar 
as  its  effect  on  the  output  of  the  array  is  concerned,  by  the  direction  of  propagation 
( e ,  4>)  of  the  wave  and  its  wavelength  X  .  The  direction  perpendicular  to  the  in¬ 
coming  wavefront  is  defined  by  direction  cosines,  cos  a,  cos/3,  and  cos  7,  where 

cos  a=  eos$  sin  6 

cos  d  =  cos$  cos  9  (1) 


cos  7  =  -sin$. 


Z 


Figure  1.  Coordinate  System  for  Arrr.y 
Structure  and  Beam  Pattern  Function 
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A  reference  plane  for  measuring  relative  phase  at  any  element  of  the  array 
[  is  chosen  to  be  that  plane  passing  through  the  origin  which  is  parallel  to  all  wave- 

fronts,  i.  e. , 

x  cos  a  +  y  cos  0  +  z  cos  7=0. 

t#  The  distance  from  any  element  to  the  reference  plane  is  then 

dj,  "  xk  cos  a  +  ^k  cos  P  4  zk  cos  y  (k  =  1,  2,  . . . ,  N)  (2) 


ana  the  phase  of  the  plane  wave  excitation  at  the  same  element  relative  to  that 
corresponding  to  the  reference  plane  is 


2t 

X 


(3) 


The  "voltage"  output  of  the  unstee.-ed  array  for  plane  wave  excitation  is  then 
given  by 


V<0,*> 


i(2*/\)dj 

L  V 

j=i 


(4) 


where  Rj  is  the  amplitude  response  ox  the  j-th  element.  Phase  delays  corre¬ 
sponding  to  preassigned  array  steering  angles  can  be  readily  incorporated  into 
the  voltage  function.  If  (0Q,  $Q)  are  the  prescribed  steering  directions,  the 
equivalent  "delay  length"  is  given  by 

1)^  =  x^  cos  aQ  +  y^  cos  /S0  *  cos  70  ,  (5) 


where  cos  «Q,  cos  cos  7Q  are  the  direction  cosines  of  the  vector  corre¬ 
sponding  to  the  steering  directions.  The  voltage  output  of  the  steered  array  is 
then 


V<0,$) 


i(2r/X)  (d:  -  DA 
Rj  e  1  J 


(6) 


“V 
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The  beam  pattern  function  in  decibels  is  defined  as  minus  the  logarithm  of 
the  ratio  of  beam  output  power  relative  to  the  maximum  output  power  over  all 
directions, 


b(0,$)  =  -10  log1Q 


V2(0,4)| 


V2(0,4)| 


0  *:  0  <  2 


T  .  T 

- <  A  <  — 

2-  -2 


(7) 


max 


In  the  event  that  the  array  is  steered,  the  denominator  in  the  previous  ex¬ 
pression  corresponds  to  the  array  response  in  the  direction  of  the  steering  vector. 

ARRAY  RESPONSE  TO  SPHERICAL 
WAVE  EXCITATION 


At  times  it  may  be  of  interest  to  examine  the  array  beam  pattern  function 
when  tie  impinging  wave  does  not  have  plane  fronts.  In  this  event,  an  alterna¬ 
tive  form  is  chosen  as  a  spherical  wave.  Thus  it  is  possible  to  examine  the  array 
responses  unde  *  other  thanfarfield  source  conditions.  Let  the  symbol  rQ  specify 
the  distance  of  a  point  source  from  the  array  center;  then  if  (0,4)  prescribes 
the  angular  orientation  of  the  source,  its  Cartesian  position  coordinates  are 


x^  =  rQ  cos4  sin  0 

y^  =  rQ  cos  4  cos0  (8) 

z^  =  rQ  sin4  . 


The  spherical  wavefront  passing  through  the  origin  is  given  by  the  equation 

(x  -  x/  +  (y  -  y^  )2  +  (z  -  z^)2  =  r2  .  (9) 

The  distance  from  an  arbitrary  element  along  a  radial  direction  from  the 
source  to  this  reference  wavefront  is  easily  found  to  be 


4 


=  r°'[ 


(xk-x£)2 1  Cvk-y{)2 


(zk  '  V 


2j1/2 


(10) 


The  ’’delay  lengths”  for  steering  directions  (0O,4O)  are  found  from  equa¬ 
tion  (5),  since  it  may  usually  be  assumed  that  the  array  beam  would  be  steered 
for  plane  wave  excitation. 
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DIRECTIVITY  INDEX 


The  directivity  index  (DI)  is  expressed  by  the  formula 

i  J-t/  2 


^2r^r/2  | 

J  J  |v2(0,*)j  cos$d$d0 


DI  =  -10  log1Q- 


REVERBERATION  INDEX 


4t  •  |v2(*,*)| 


max 


The  reverberation  index  (RI)  is  expressed  by  the  formula 


RI  =  - 10  log 


nRV(M)  |v2(0,«)|  cos*d<i>d0 

r/9.  1  1 


10  ‘ 


4r*  |v2(0,$)| 


max 


(11) 


(12) 


where  RV(0,$)  is  a  reverberation  function  which  is  used  to  effectively  describe 
the  relative  reverberation  strength  of  backscattered  energy  from  boundary  or 
volume  scatterers  over  the  angular  intervals  involved. 


COMPUTATIONAL  ASPECTS 

The  computer  program  developed  to  perform  array  beam  calculations  con¬ 
tains  three  separate  sections: 

1.  Determination  of  the  array  beam  pattern 

2.  Evaluation  of  DI 

3.  Evaluation  of  RI  for  an  arbitrary  reverberation  function. 

Given  ihe  physical  structure  of  the  array,  the  first  of  these  merely  involves  per¬ 
forming  the  summation  given  in  equation  (6)  for  either  the  plane  or  the  spherical 
wave.  Provision  is  made  to  insert  the  individual  element  response  and  a  possible 
coefficient  for  amplitude  shading  from  two  separate  subroutines  of  the  CALL 
statement.  Similarly,  the  physical  delays  as  well  as  those  artificially  introduced 
for  array  steering  are  computed  in  another  such  subroutine.  The  computed  beam 
function  is  presented  as  either  a  polar  or  a  rectangular  plot,  as  required  by  the 
user.  Details  of  input  information  and  use  of  the  program  are  given  in  later 
sections. 
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The  second  function,  computation  of  directivity  index,  involves  integration 
of  a  two-dimensional  integral  of  the  form 


1  = 


f(x,y)  dxdy  . 


(13) 


The  numerical  evaluation  of  the  above  integral  is  accomplished  by  means  of 
a  quadrature  formula. 


M  N 

Iss£  S  HiHj  f«*i.yj)  •  (i4) 

i=i  j=i 

The  thi  rd  program  fimction  is  similar  to  the  above,  since  computation  of  the 
reverberation  index  also  involves  the  quadrature  formula.  In  addition,  it  is  nec¬ 
essary  to  provide  a  reverberation  function  which  is  to  be  generated  in  a  subroutine 
provided  by  the  user. 

The  integration  mesh  size  is  a  two-dimensional  variable  in  order  to  econo¬ 
mize  on  computation  time.  In  the  absence  of  round-off  errors,  a  sufficiently 
small  mesh  size  enables  the  quadrature  formula  to  adhere  to  a  certain  desired 
computation  accuracy. 

DEFINITION  OF  INPUT  PARAMETERS 


The  input  parameters  appearing  on  the  NAMELIST  of  the  executive  program 
have  four  different  categories: 

1.  ARRAY  CHARACTERISTICS 

2.  BEAM  PATTERN  FUNCTIONAL  CONTROL  PARAMETERS 

3.  DI  AND  RI  COMPUTATIONAL  CONTROL  PARAMETERS 

4.  CONTROL  PARAMETER  OF  SEQUENTIAL  RUNS. 

They  are  defined  below. 

1.  ARRAY  CHARACTERISTICS 
DIMENSION  INDEX 

L  Number  of  elements  in  the  x-curectioii 

M  Number  of  elements  in  the  y-direction 

N  Number  of  elements  in  the  z-direction 

LMN  Total  number  of  elements  (for  irregular  element 

spacing) 
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FREQUENCY  (Hz) 

F  \  Frequency 

SOUND  VELOCITY  C 

VE1|CTY  If  unspecified,  the  ’’default"  value  is  4900.0  ft/sec 
SPACING  E-ETWEEN  ELEMENTS 


LXYZ  =  0  Implies  elements  equally  spaced  in  x,  y,  and  z 
directions 

XSPACE  Spacing  between  elements  in  x-direction  (ft) 

YSPACE  Spacing  between  elements  in  y-direction  (ft) 

ZSPACE  Spacing  between  elements  in  z-di  recti  on  (ft) 

♦♦A  nonzero  indicator  IXYZ  requires  that  each 
element  location  be  given  in  three-dimensional 
Cartesian  coordinates.  These  can  be  given  in  free 
format. 


SELECTABLE  PHASE  DELAY  OPTION 

ID  ELAY  =  0  Implies  no  preselected  phase  delay  associated  with 
elements. 

♦♦If  IDELAY  is  not  0,  DELAY  subroutine  should  be 
supplied  with  selectable  delay  information 

ARRAY  ELEMENT  DIRECTIONALITY 

iR  =  0  Implies  nondirectional  elements;  otherwise  the  sub¬ 

routine  RSPNSE  has  to  be  supplied,  containing  the 
Information  of  directionality  in  the  array  R 


SHADING  TECHNIQUE 


IS  HAD  =  0  Implies  no  shading  cf  individual  element. 

♦’•'Nonzero  IS  HAD  requires  an  amplitude  shading 
technique  to  be  given  in  SHADE  subroutine 


WAVEFRONT 

RZERO  <  0  Plane  wavefront 
RZERO  >  0  Spherical  wavefront 

♦♦RZERO  is  the  distance  from  a  source  to  the  origin 
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MAXIMUM  ARRAY  RESPONSE  ON  MAJOR  AXIS 

MRA  Set  =  0  initially. 

The  steering  direction  must  be  given  to  determine 
the  maximum  response.  Once  the  maximum  response 
value  is  determined,  the  program  automatically  sets 
MRA  / 0.  For  subsequent  runs,  the  latest  maximum 
response  value  is  used,  unless  a  different  steering 
direction  is  specified  with  MRA  =  0. 

THEM  Horizontal  steering  direction  angle 

PHIM  Vertical  steering  direction  angle 

2.  BEAM  PATTERN  FUNCTION  CONTROL  PARAMETERS 

BEAM  TYPE 


IBEAM  =  0 
IBEAM  =  1 
IBEAM  >  1 

PLOTS 

IPLOTcO 
IPLOT  =  0 
IPLOT  >0 
PLOT  SCALES 
ALONG 
DLLONG 
DIAM 
PLDB 


Horizontal  beam 
Vertical  beam 

Calculates  major  response  value;  implies  no  need 
to  calculate  the  beam  pattern;  calculations  proceed 
directly  to  the  computation  cf  DI/RI 


Plot  is  suppressed 

Rectangular  plot 

Both  rectangular  and  polar  plots 

Rectangular  [dot,  horizontal  length  in  inches 
Rectangular  plot,  vertical  scale  in  decibels 
Polar  plot,  radius  in  inches 
Polar  plot,  scale  in  decibels 


CONTROL  OF  NUMERICAL  OUTPUTS 
OF  BEAM  PATTERN  FUNCTION 


IPRINT  =  0  Printout  of  numerical  values  is  suppressed 
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ANGULAR  LIMITS  AND  INCREMENTS  (Degrees) 


THMIN 

THMAX 

PHIMIN 

PHIMAX 

DTHETA 

DPHI 


Minimum  horizontal  angle  for  beam  computation 

Maximum  horizontal  angle  for  beam  computation 

Minimum  vertical  angle  for  beam  computation 

Maximum  -  rtical  angle  for  beam  computation 

Horizontal  angular  increment  between  successive 
computations 

Vertical  angular  increment  between  suecessi  «/e 
computations 


3.  DI  AND  RI:  COMPUTATIONAL  CONTROL  PARAMETERS 
CONTROL  OF  DI  COMPUTATION 


INDEX  ~  0  Implies  that  the  calculation  of  DI  is  to  be  omitted 


CONTROL  OF  RI  COMPUTATION 

IREV  t  0  Hie  calculation  of  I J  is  required.  The  reverberation 

function  RV  ($,$)  must  be  given  in  .•»  subroutine  REVB 


INTEGRATION  LIMITS  AND  SPACING 

HHZ  Horizontal  integration  step  size 

Default  value  =  3° 

HVT  Vertical  integration  step  size 

Default  value  =  3° 

LOWH  Lower  limit  of  integration  in  horizontal  direction 

Default  value  =  0 

LIMH  Number  of  points  considered  between  horizontal 

limits 

Default  value  =  120 

LOWV  Lower  limit  of  integration  in  vertical  direction 

Default  value  =  -90* 

Number  of  points  considered  between  vertical  limits 
Default  value  -  60 


LEMV 
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INTEGRATION  CONTROL 


i: 


MANYH  Number  of  intervals  to  be  integrated  using  different 
step  sizes 

Default  value  =  1 

4.  CONTROL  PARAMETER  OF  SEQUENTIAL  RUNS 

ITIME  =  0  The  very  first  time  entering  the  program.  (The  pro¬ 
gram  initializes  th?  quantity  to  be  0) 

ITIM2  =  1  Implies  there  is  no  need  to  repeat  the  calculations 
of  the  array  elements,  the  element  response,  delay, 
shading,  and  the  Major  Response. 

ITIME  >  1  Implies  there  is  no  need  to  calculate  the  beam  pattern 
function,  plots,  and  related  computations. 

Other  than  the  control  parameter  of  sequential  runs,  the  foregoing  tabulation 
of  input  parameters  is  summarized  in  tables  1,  2,  and  3. 


Table  1.  Array  Characteristics 


n^criptica 

Dimension 
Number  of  element* 

Frequency 

Velocity 


j  Element*  (pacing 


Element*  response 


Shading  technique 


Wavefront 


Symbol 

L,  M,  N 
IMS 


VELCTY 

IXYZ 

XSPACE 

YSPACE 

ZSPACE 

IDELAY 

DELAY 

DELA 

1R 

RSPNSE 

R 

IS  HAD 
SHADE 
SHAD 

R/.ERO 


Default  Condition 


<K)0.0  ft/ sec 


Array  response  on  major  axis 
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PROGRAM  DESCRIPTION 

The  overall  structure  of  the  computer  program  is  described  by  means  of 
table  4,  the  following  brief  outline  of  the  executive  program,  and  the  various 
subroutines  controlled  by  the  program. 

Table  4.  Structure  of  the 
Computer  Program 

GBEAM 

- SUM 

- XYZ 

- RSPNSE 

- DELAY 

- SHADE 

- BPLOTR 

- BPLOTP 

I - RADIAL 

i - REVE 

GBEAM  —  This  is  the  executive  program  that  commands  the  specific  perform¬ 
ance  requested  by  the  user.  It  performs  the  following  functions: 

1.  Accepts  input  data 

2.  Calculates  the  wavelength 

3.  Calls  RSPNSE,  DELAY,  SHADE,  if  they  are  required 

4.  Commands  the  setup  of  clement  spacing,  if  equal  spacing  is  required 

5.  Lists  certain  input  information 

6.  Displays  some  processed  information 

7.  Determines  Major  Response  of  array 

8.  Calls  for  rectangular  and/or  polar  plots 

9.  Calculates  horizontal  and/or  vertical  beams 

10.  Calculates  DI 

11.  Calculates  RI 

12.  Prints  out  numerical  results. 
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The  various  subroutines  are: 

SUM  Calculates  Major  Response  and  the  delays 

XYZ  Calculates  the  element  locations,  if  equally  spaced,  and  accepts 
irregularly  spaced  array  element  locations  provided  as  input  to 
the  program 

RSPNSE  Obtains  element  response,  if  elements  do  not  possess  an  omui di¬ 
rectional  character 

DELAY  Accepts  selectable  delays 

SHADE  Accepts  an  arbitrary  amplitude  shading  technique 

BP  LOIR  Generates  a  rectangular  plot 

BPLOTP  Generates  a  polar  plot 

RADIAL  Sets  up  radial  scale  for  polar  plot 

REVB  Describes  a  reverberation  function. 

PROGRAM  USE 

The  coding  of  the  program  for  the  Univac  1108  computer  is  compatible  with 
the  EXEC  8  system.  The  program  has  >  een  labeled  GBEAM  and  is  availabl  e  on 
tape  (Z158,  1st  file).  Under  ordinary  use,  no  additional  tapes  are  needed.  The 
input  parameters  previously  defined  should  be  specified  in  NAMELIST  format 
under  the  name  INPUTS,  and  the  input  procedure  should  therefore  present  no 
great  difficulty.  Several  points,  however,  require  explanation. 

to  the  event  that  the  array  spacing  is  irregular,  the  proper  indication  must 
be  given  (IXYZ  >0)  in  Hie  NAMELIST  data.  When  all  the  NAMELIST  data  have 
been  entered,  the  specific  locations  of  all  array  elements  in  Cartesian  coordi¬ 
nates  with  the  ordering  x,  y,  z  should  be  entered  in  free  format.  When  the 
executive  program  GBEAM  calls  for  a  specific  element  position  information, 
the  subroutine  'TZ  reads  the  required  information  and  passes  this  on  as  re¬ 
quired.  (When  <. —  -ay  elements  are  regularly  spaced,  the  subroutine  XYZ  calcu¬ 
lates  file  coordinates  of  the  element  specified  from  the  given  array  spacing 
information. ) 

The  program  contains  a  number  of  other  optional  features  which  require 
special  treatment  by  the  user.  Individual  elenr  sots  may  be  given  an  amplitude 
response  other  than  omnidirectional  by  use  ot  the  subroutine  RSPNSE,  in  which 
the  user  must  define  (for  example,  by  means  of  an  equation  or  tabular  data)  the 
desired  response  function  as  it  depends  upon  the  deviation  angles  9  and 
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/ 


Similarly,  an  ampi’tude  coefficient  compatible  with  some  desired  shading 
technique  to  the  beam  may  be  introduced  as  a  modification  to  the  individual  ele¬ 
ment  response  by  the  subroutine  SHADE.  An  example  which  serves  to  illustrate 
the  manner  in  which  SHADE  subroutine  may  be  used  is  given  in  the  appendix. 


The  subroutine  DELAY  calculates  the  artificial  delay  lengths  associated 
with  a  given  tilt  and  steering  angle.  Also,  it  is  possible  to  introduce  an  arbi¬ 
trary  phase  delay  into  the  response  of  each  element  by  providing,  as  input,  the 
required  numerical  values  in  free  format. 


The  basic  options  of  die  program  are  shown  in  table  5,  which  defines  the 
various  input  quantities  according  to  the  option  desired.  If  a  reverberation  index 
is  desired,  the  subroutine  REVB  must  be  supplied  with  the  reverberation  function 
RV(0,<*>)  in  a  manner  analogous  to  that  used  in  RSPNSE,  which  was  previously 
discussed. 


Table  5.  Control  Parameters 


Geaenl  Optic 


•  Items  required 
i.  Items  to  be  specified 
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OUTPUT  INFORMATION 

The  lirst  page  of  the  output  is  given  as  an  aid  to  the  user  ir.  cheeking-  the 
input  data  and  consists  of  the  following  information: 

Title:  GENERALIZED  BEAM  PATTERN  PROGRAM 
Wavefront  type 
Array  size 

Frequency  (Hz)/wavelength 
Spacing  in  all  directions  (ft) 

Reference  point  (geometric  center  of  the  array) 

Selectable  phase  delays  introduced 
Deter  min  id  array  response  on  major  axis. 

Subsequent  pages  contain  horizontal  or  vertical  beam  pattern  numerical  out¬ 
puts  if  the  printout  optical  is  selected. 

The  next  output  information  gives  the  numerical  values  of  Dl/RI  if  requested. 

The  last  item  of  output  information  indicates  how  many  plots  will  be  generated 
if,  indeed,  plots  have  been  required. 

IMPORTANT  FEATURES 
Simple  mathematical  formulas. 

Program  is  generai,  complete,  optional,  easy  to  use,  and  has  no 
restrictions. 

Important  parameters  are  well  defined  internally.  No  danger  is  involved 
if  the  user  failed  to  specify  some  of  them.  In  addition,  the  user  has  the 
optical  of  altering  these  parameters  to  satisfy  his  requirements. 

Capability  of  handling  both  spherical  and  plane  wavefronts. 

Rectangular  and  polar  plots  are  built  in.  User  can  select  the  plot  size. 

Automatic  setup  of  equal  array  element  spacing. 

Accepts  irregularly  spaced  array  elements  by  their  locations. 
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The  use  of  steering  direction  information  to  determine  the  MRA 
results  in  considerable  time  saving. 

Option  of  accepting  the  element  directionality. 

The  SHADE  subroutine  accepts  any  shading  technique. 

The  program  arranges  the  three-dimensional  elements  in  one¬ 
dimensional  computations. 

Integration  limits  are  adjustable  so  that  the  computations  can  be 
accelerated.  Program  can  decompose  a  two-dimensional  integral 
into  a  set  of  subintegrals  along  each  subintegral  with  different  step 
sizes;  the  step  size  can  be  any  fraction  of  one  degree. 

The  listing  of  the  basic  input  information  on  the  first  output  sheet 
enables  the  user  to  examine  whether  or  not  his  input  information  is 
supplied  correctly. 

All  subroutines  need  not  be  supplied  by  the  user.  If  they  are  not  called 
for,  the  program  takes  care  of  them  automatically. 

Manipulation  cf  the  large  number  of  array  elements  is  flexible  to  the 
upper  limit  of  computer  capacity. 

The  option  of  computing  and  plotting  any  one  section  of  the  beam 
pattern. 

The  option  of  bypassing  the  beam  plots  to  calculate  the  DI  and/or  RI. 

If  any  additional  computations  are  required,  the  necessary  modifica¬ 
tions  can  be  incorporated  easily. 


NUMERICAL  EXAMPLES 


In  this  section,  selected  examples  are  presented  to  demonstrate  the  options 
of  the  program  and  to  illustrate  typical  computation  results.  To  demonstrate 
subroutine  usage,  a  cos^  shading  technique  is  applied  in  one  problem.  The 
examples  are  divided  into  three  categories.  The  three  arrays  considered  have 
the  configurations  described  in  figure  2. 
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Categor>r  1  —  the  first  division  has  four  examples  of  one-dimensional,  equally 
spaced  elements  (table  6).  Initially,  a  plane  wavefront  is  consid¬ 
ered.  Following  this,  three  spherical  wavefronts  are  presented: 
one  has  a  short  distance  from  the  source  to  the  array  center;  one 
has  a  moderate  source  distance;  and  one  has  a  great  source  dis¬ 
tance.  The  effects  of  the  three  different  source-array  separation 
distances  are  shown  graphical  >y  in  figures  3  through  6.  Each  fig¬ 
ure  consists  of  rectangular  and  polar  plots  of  the  horizontal  and 
vertical  beams  for  each  example.  The  plots  in  the  figures  are  iden¬ 
tified  numerically  with  their  corresponding  example  numbers  in 
the  table. 

Category  2  —  A  two-dimensional,  equally  spaced  rectangular  array  with  a  plane 
wavefront  is  considered  next  (table  7).  Initially,  the  beam  is  com¬ 
puted  without  shading,  then  with  a  cos2  shading  applied  in  two 
dimensions.  Each  example  in  the  table  has  two  different  plots. 
Figures  7  and  8  each  contain  four  plots  of  either  the  horizontal  or 
vertical  beam  p^^eras  for  the  rectangular  and  polar  plots  of  each 
example. 


Table  6.  Category  1  —  One- Dimensional  Array  (9x1x1) 


Example  Number 

1.1 

1.2 

1.3 

1.4 

Array 

Excitation 

Type 

Plane 

Wave 

Spherical 

Wave 

Spherical 

Wave 

Spherical 

Wave 

Source 
Distance  (it) 

Far 

Field 

5.0 

20.0 

1000. 0 

1 

Element  Spacing  (X) 

0.5 

0.5 

0.5 

0.5 

Horizontal 

Beam 

R-plot* 

plot  1. 1. 1 

plot  1.4.1 

P-plotf 

plot  1. 1.3 

plot  1.2.3 

plot  1.3.3 

U 

A 

R-plot 

plot  1. 1. 2 

plot  1.2.2 

plot  1.3.2 

plot  1.4.2 

P-plot 

plot  1.1.4 

plot  1.4.4 

Directivity  Index 

9.54 

9. 65 

9. 55 

9.54 

Reverberation  Index 

12.55 

12. 66 

12. 56 

12.55 

Computation  Time  (sec) 

27 

30 

31 

32 

•R-plot  ^  rectangular  plot 
tp-plot  -  polar  plot 
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Table  7.  Category  2  —  Two-Dimensional  Array  (20  x  1  x  3) 


Example  Number 

2.1 

2.2 

Shading  Technique 

None 

coe2 

Array 

Excitation 

Type 

Plane  Wave 

Plane  Wave 

Source 
Diatance  (ft) 

Farfleld 

Farfleld 

Element  Spacing  (X) 

0.5 

0.5 

Horizontal 

Beam 

R-plot* 

plot  2. 1. 1 

plot  2. 2.1 

P-plott 

plot  2. 1.3 

plot  2. 2. 3 

Vertical 

Beam 

R-plot 

plot  2. 1.2 

plot  2.2.2 

P-plot 

plot  2. 1.4 

plot  2.2.4 

Directivity 

Index 

Fixed 

Step  Size 

19.03 

Variable 

Step  Size 

19.08 

19.00 

Reverberation 

Index 

Fixed 

Step  Size 

22.04 

Variable 

Step  Size 

22.09 

22.01 

Computation 
Time  (tec) 

Fixed 

Step  Size 

101 

Variable 

Step  Size 

58 

99 

'Rectangular  plot 
t Polar  plot 
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Figure  8.  Vertical  Beam  Patterns  of  Rectangular  Array 
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Category  3  —  Finally,  a  three-dimensional  cylindrical  array  of  21  elements  with 
different  wavefronts  is  considered  (table  8).  The  cylindrical 
height  =  IX,  the  vertical  spacing  0. 5X,  and  the  circumferential 
spacing  is  (ir/6)X,  or  roughly  0. 5X,  and  the  radius  is  IX.  This 
category  demonstrates  the  treatment  for  irregularly  spaced  ele¬ 
ments.  In  table  8,  each  example  has  two  different  plots.  The  four 
plots  of  figures  9  and  10  show  either  tho  horizontal  or  vertical 
beam  patterns  for  the  rectangular  and  polar  plots  of  each  example. 


In  the  examples,  whenever  an  RI  calculation  is  involved,  the  function  RV(0,4) 
is  defined  for  illustrative  purposes  to  be  a  constant  0. 5  for  all  0  and  4.  There¬ 
fore,  no  physical  significance  is  intended.  Steering  directions  0  =  4  =  0  are  used 
for  examples  of  categories  1  and  2;  0  =  45*,  4=0  are  used  for  examples  of 
category  3.  The  distance  between  elements  is  chosen  to  be  half-wavelength,  if 
equal  spacing  is  required.  Computations  of  DI  and  RI  are  given  in  tabular  form, 
together  with  the  computation  time. 


Table  8.  Category  3  —  Three-Dimensional  Array  (21  Elements) 


•  Example  Somber 

3.1 

3.2 

Array 

Excitation 

Type 

Plane  Wave 

Spherical  Wave 

Source 

Distance  (ft)  i 

Farfield 

5.0 

Element  Spacing 

Irregular 

Irregular 

Horizontal 

Beam 

R-plot* 

•at  3. 1.  1 

plot  3.2.1 

P-plott 

plot  3.1.3 

plot  3. 2. 3 

Vertical 

Beam 

R-plot 

plot  3. 1. 2 

plot  3.2.2 

P-pIot 

plot  3. 1.4 

plot  3.2.4 

Directivity  Index 

13.02 

12.80 

Reverberation  Index 

16.03 

15. 81 

Com  potation  Time  (aec) 

57 

54 

•Rectangular  plot 
t Polar  plot 

26 


01  DOWN 


TR  4687 


HCl  3.1.3  Cnir.D«'CM  aiut.  n*f. I  w*vlMO.: 

U mV  »A1!i»N 

•  0.0 


not  3.JJ  ttCtANCUA*  A*«AV  V^I*«CA|  WAVtltC^T 
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9  0.0 


Figure  9.  Horizontal  Beam  Patterns  of  Cylindrical  Array 
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Figure  10.  Vertical  Beam  Patterns  of  Cylindrical  Array 
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Example  2. 1  of  table  7  has  two  values  of  DI  and  HI,  which  are  obtained  by 
using  two  different  mesh  sizes  for  the  numerical  integration.  The  value  above 
the  line  is  calculated  by  integrating  horizontally  over  the  interval  (0°,  90°),  using 
1°  spacing;  the  value  below  the  line  is  calculated  by  the  summing  of  two  integrals. 
One  is  integrated  over  the  interval  (0®,  15*)  with  a  0. 5°  mesh  size;  the  other  is 
integrated  over  the  interval  (15*,  90*)  with  a  1. 5®  spacing. 

REMARKS 

The  present  limit  on  the  total  number  of  elements  is  2400.  The  number  of 
the  integration  points  must  be  an  odd  number,  i.  e. ,  the  LIMH,  LIMV  must  be 
an  even  number  >2. 


REFERENCES 

1.  P.  Davis  and  P.  Rabinowitz,  Numerical  Integration,  Blaisdell  Publishing  Co. , 
New  York,  1967. 
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3.  R.  J.  Urick,  Principles  of  Underwater  Sound  for  Engineers,  McGraw-Hill 
Bode  Co. ,  Inc. ,  New  York,  1967. 
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Appendix 

COMPUTER  PROGRAMS 


This  section  lists  all  FORTRAN  V  computer  programs  required  for  the 
computations. 

Subroutines  KSPNSE  and  DELAY  are  made  dummies  to  satisfy  the  compiler 
for  testing  purposes. 

Subroutine  SHADE  describes  a  pedestal  shading  technique  that  is  demonstrated 
in  one  of  the  test  examples. 
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GBEAM 

C  «***  A  GENERALIZED  BEAM  PATTERN  COMPUTATION  PROGRAM 

C  WITH  SHADING*  directivity  index  and  reverberation  level 

C  ****  D.lEE  FEB.  1973  -  SYSTEMS  ANALYSIS  DEPARTMENT 
c  ****  IRREGULAR  LOCATIONS  of  ELEMENTS  IN  3-D I MENS I ON 
C  ♦*♦♦  IR=0  OMNI-DIRECTIONAL 
C  ♦•♦*  IXyZ=0  ELEMENTS  EQUALLY  SPACED 
C  ♦♦♦♦  FOR  SPHCRiCAL  WAVE  FRONT*  ENTER  NONZERO  RZERO 
C  *»**  IF  NOT  REUUESTLD,  HORIZONTAL  BEAM  IS  GIVEN 
C  IBEAM=1  SVERTICAL  BEAM 

DIMENSION  R(2400) ,XH(24Q0) » YH12400) ,ZH(2400) 

DIMENSION  DEL A (2400) » SHAD (2400 ) » BEAM ( 720 ) 

DIMENSION  PANGLEC 720) .TANQLfc (720) 

DIMENSION  HUR(12)*DATA(200) 

CALL  MODESG(DATA.O) 

1  FORMAT (iOX»3(E15»8.15X) ) 

2  FORMAT (SIX. 'GENERALIZED  BEAM  PATTERN  PROGRAM'//) 

3  FORMAT (/bX. 'ARRAY  SIZE',l<iX,'(«,l3,lX,'X*,lX,X3,lX,'Xt,lX,I3»*)*) 

4  FORMAT (/5X»  'FREQUENCY*  »UX#  's' ,El5»8,2X,  •  (HZ)  • ) 

5  FORMAT (/bX, 'SPACING  Z-OlRgCTION  s',tl5.8»2X* ' (FEET) • ) 

6  FORMAT (/SX» 'SPACING  X-OXRECTION  =»»E15,8r2X* * (FEET) • ) 

7  FORMAT (/5X, *****  HORIZONTAL  BEAM  ♦♦♦♦•//) 

B  FORMAT (/SX. *****  VERTICAL  BEAM  ****«//) 

9  FORMAT (/bX,* SELECTABLE  PHASE-DELAY  IS  ABSENT') 

10  FORMAT (1H1) 

11  FORMAT(/5X» 'SPACING  Y.OIRECTION  =» »EX5.8*2X, • (F£ET) ' ) 

12  FORMAT (/5X. 'SPACING  IN  X-OlRECTlON  IRREGULAR') 

13  FORMAT (/SX? 'SPACING  IN  Y-DIRECTION  IRREGULAR*) 

14  FORMAT (/bX, 'SPACING  IN  Z-OXRECTION  IRREGULAR' ) 
lb  FORMAT (40X,2(E15.8,15X) ) 

17  FORMAT (/SXf'MRA  GIVEN  AT  (TmETA”'»F6,1,'  DEGREES*  PHl=»,F6.1»'  DEG 
♦REES) ' ) 

18  FORMAT (12A6) 

19  FORMAT (/SX* 'PLANE  WAVE  FRONT') 

20  FORMAT (/SX* 'DIRECTIVITY  XNOEX  =»,F8.2,2X» ' (DB) » ) 

21  FORMAT (//SX, 'REVERBERATION  LEVEL  s»,F8,2»2X* » (DB) • ) 

501  FORMAT (/5X, 'SPHERICAL  YAVE  FRONT  —  R0=* »El6.8,2X, • (FEET) • ) 

502  FORMAT (/5X» 'LOCATIONS  CF»*l4,lX»'AKKAY  ELEMENTS  ARE  GIVEN') 

503  FORMAT (/5X» 'REFERENCE  POINT  ( » ,F6,2,1X*' » » ,F6.2,iX» ' >  •  »F6*2,  •)  •) 

504  FORMA T(/bX» 'WAVE  LFJ^GTH  s» »F6.2,3X, » (FT) ' ) 

NAMELIST  inputs/ir, ixyz»l*m,n» xdelat *f» rzero, 

♦  XSPACE  * YSPACE* ZSPACE* DTHETA  *  DPMI *  I BEAM »  THEM , 

♦  phim»velcty,phimin»phxmax,thmin,thmax, 

♦  I PLOT, INDEX* IREV* ITIME* ALONG, DLLONG*ISHaO, 

♦  H8Z*HVT'XPRINT*L0WH*LXMh*L0WV,LIMV'MAnYH/ 
data  IR* IXYZ, IPcOT, XBEAM* INDEX, IREV* ITXMEiMRA, ISHA0/9*0/ 

DATA  MANYH*lPRlNT,LOWH»LlNH,LOWV*LiMV/l,2wO* 120,-90*60/ 

DATA  HHZ*HVT/2w3»0/ 

OATA  PHlMlN»PHIMAX»rHMXN,THMAX»VELCTY/-90, 0*90, 0*0, 0*360*0*4900./ 
OAT A  OTHCTA  r  OPHI *  PI , OEG/2* 10.0*3,1415926535*. 01745329/ 

DATA  ALONG, D-LONG,PLDB»DXAM/2#6,0*60. 0*4.0/ 

OAT A  R» SHAD, DEL A/480 Owl , 0*240040 • 0/ 

100  READ  (3, INPUTS, ERR=9999,EN0s9999) 
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REAO  18.  HUH 

IFdTlME  ,NE.  0)  60  TO  86 
WLAmDA=V£UCTY/F 

CAUL  XYZUXYZ.XH.YH.ZH.L.M.N.XSPACE.YSPACE.ZSPACE.LMN.XO.YO.ZO) 
IFdR  .NE.  0}  CAUL  RSPN5E(LMN.R> 

1F(10ELAY  ,NE»  0)  CAUL  DELAY ( UMN . DELA ) 

IFdSHAO  »Nt,  0}  CALL  SHADE(LMN.XH» YH.ZH.XO.YO.ZO) 

RMAX=0,0 
00  89  Isl.LMN 

1FCRMAX  «LT •  R(l))  RMAXsR(I) 

89  CONTINUE 
88  CONTINUE 

IFdTIHE  ,6T,  I)  GO  TO  300 
IFdTIKE  ,£tt,  1)  60  TO  1111 
WRITE  (4.2) 

IF(RZERO  .Lt.  0.0)  WRITE  (6.19) 

IFCrZERO  .»T.  0*0)  WRITE  (4.501)  RZERO 
IF(l  »N£.  0)  WRITE  (4,3)  (..m.N 
IF(l  .EQ.  0)  WRITE  (4.502}  lMN 
IF (F  .LE.  0.0)  WRITE  (4.SQ4)  WUAMQA 
IF(F  .6T.  0.0)  WRITE  (4.4)  F 

IFdXYZ.EQ.O  ,ANO.  XSPACE.6T.0.0)  WRITE  (4.6)  XSPACt 
lFdXYZ.ttt.O  ,ANO,  YSPACE.6T.0.0)  WRITE  (4.H)  YSPACE 
IFdXYZ.EQ.O  .AND.  ZSPACE.6T.0.0)  WRITE  (4.5)  ZSPACE 
WRITE  (4.503)  XO.YO.ZO 
IFdOEUAY  ,EQ.  0)  WRITE  («»9) 

IF(MRA  .NE.  0)  60  TO  52 

CAUL  SUM(Th£M.PHIM»RZERO»XH.YH.ZH»WLAHDA'LMN.DElA.R.RmAX»ARm.MRA»S 
•HAD) 

52  WRITE  (4.17)  THEM.PHIM 
WRITE  (4.10) 

1111  continue 

IFdBEAM.EU.O  .AND.  IPRINT.EQ.O)  WRITE  (4.7) 

IFdSEAN.EQ.l  .AND.  IPRINT.EQ.O)  WRITE  (4.8) 

IFdBEAH  «EQ.  0)  60  TO  200 
IFdBEAM  ,6T.  1)  60  TO  300 
C  *••*  VERTICAL  BEAN 

theta=thmjn 

JNDEX=0 

30  JN0EX=JNDEX^1 

PHIsPHIMIN+FLOAT ( JNOEX-1 ) WOPHI 
PAN6UE(JN0EX)SPHI 
IF(PHI-PHIHAX)  32.32.38 
C  *•*•  FOR  A  FlXtO  THETA.  PUOT  PHI  VS  BEAM 
38  JNP£X=JN0EA-1 

IF(IPU0T  ,LT.  0)  60  TO  300 

CAUL  BPLOTK ( IBEAM * PAN6UE .BEAM. JNOEX » THETA » PHIMIn . PHIMaX . ALONG . DUUO 
•N6.mDR.QATA) 

IFdPLOT  »ST»  0)  CAUL  BPLOTP(IBEAM>PANGLE»BEAM»jNDEx»THETA»PLDB»DI 
•AM.h0R.DATA) 

60  TO  300 

32  CAUl  SUM (THtTA.PHI .RZERO. XH,YH»ZH»WLAMOA.UMN»DElA»R.RmAX»XMAX»MRA. 
•SHAu) 

beam ( JNOEX ) =-10 . • AL0410 ( XMAX/ARM ) 

IFdPRINT  .NE,  0)  60  TO  30 
IF (JNOEX  .EG.  1)  60  TO  70 
WRITE  (4.15)  PHI. BEAM (JNOEX) 

60  TO  30 


33 


TR  4687 


N  ~ 


70  WRITE  (4»1>  THETA » PHI, BEAM (JNOEX) 

60  TO  30 

€  •***  HORIZONTAL  BEAM 
200  CONTINUE 
PHUPHImIN 
JN0£X=0 

210  JNO£X=JNOEX+l 

THETASTHM Jn*FLOAT ( JNOEX-1 ) *OTHETA 
TAN6LE (JNOEX) =THETA 
IF (THETA-7HHAX)  213*213*214 
C  ***•  FOR  A  FIXED  PHI*  PLOT  THETA  VS  BEAM 

216  JNO£X=JNUEX-l 

IFUPLOT  ,LT.  0)  60  TO  300 

CALL  BPLOTK ( ISEAM * TAN6LE ,  BEAM , JNOEX  »PHX  , THMIN,  ThMAX  » AL0N6 .0LL0N6 , H 
*UR»OATA) 

IFUPLOT  ,6T«  0)  CALL  BPLQTP(XBEAM*TAN6lE*BEAM» jNO£X*pHl*PLOB*DlAM 
**hdr*dataj 

60  TO  300 
213  CONTINUE 

CALL  SUM(THETA*PHlfftZERO'XH»YH»ZH'HLAMDAtLMN'DELA»R'RMAX’XnAX'MRA' 
•SHAO) 

BEAM (JNOEX ) =-10 . *AL0610 ( XMAX/ARH) 

XFUPRZNT  .NE.  0)  60  TO  210 
IF( JNOEX  «C6»  1)  60  TO  230 
WRITE  (4»l&)  THETA* BEAM (JNOEX) 

60  TO  210 

230  WRITE  (4*1)  PHI, THETA, BEAM(JNOEX) 

60  TO  210 

300  1FUPRINT  .EO.  0)  WRITE  (4*10) 

C  TO  CALCULATE  DIRECTIVITY  ZNOEX  *  REVERBERATION  LEVEL 
IFUNOEX.E6.0  »ANO»  IREV.EO.O)  60  TO  100 
TO 1=0,0 
TRLsO.O 

1210  IF(MANYH  *LE.  0)  60  TO  1200 

READ  (3*INPUTS*ERR=9999*EN0s9999) 

F I HoZ=HHZ* JE6 

FIVT=HVT*0E6 

HLOw=FLOAT(LOWH) 

VLOwsFLOAT(LOWV) 

OINOEX=0*0 

RLEVEL=0.0 

1IN0=U 

00  1000  JO=0*LXMH 
01=0.0 
RL=0.0 

TH=FLOAT (JO)*HHZ+HLOW 
JINO=+i 

00  1001  KO=0*LXMV 
PH=FLOAT  (KO)  *HVUVLOW 

CALLSUM(TH»PH*RZERO*XH»YHeZH*WLAMOA*LMN,OELA*R*RMAX*XMAX»MRA*SHAD) 
ZFUREV  .NE.  0)  CALL  REVB(REV»TH'PH) 

X=XmAX*C0S(PH«D£6)/ARM 
IFUREV  .NE.  0)  YsREVwX 
XF(KO  .NE,  0)  60  TO  1002 
1003  OI=Ol+X 

IFUREV  .NE.  0)  RL=RL4Y 
60  TO  1001 

1002  IF(KO  .£6,  LIMV)  60  TO  1003 
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ox=oi+x+x 

XFUREV  .NE.  0)  RL=RL*Y+Y 
IFtJlND  .GT.  0)  l>X=OI*X+X 
IFtjXND.GT.Q  .AND,  IREV.NE.O)  RL=RL+Y+Y 
JINO=-JINO 
1001  CONTINUE 

IF (JO  ,Nt#  0)  60  TO  1005 

1004  DIN0EX=0lN0EX+FiVT*DI/3. 

XFUREV  ,NE,  0)  RLEVEL=RLEVEL+FIVT*RL/3. 

60  TO  1000 

1005  1F(J0  ,£«,  LIMH)  60  TO  1004 
0lN0£XsDlNUEX+2,*FIVT*DI/3. 

XFUREV  «N£,  0)  RLEVEL=RLEVEL+2.*FIVT*Rl/3. 

IFUXNO  «&T j  0)  0INUEX=0XN0tX+2.*FIVT*DX/3. 

xfuino.gt.o.ano.xrev.ne.o)  RLEVEL=KLEVEL+2.*FIvT*RL/3. 
XXNqs-XXNO 
10Q0  CONTINUE 

TOIsT0X4FXHO2*OINOEX/S«0 
TRLsTRL+F inO£*RLEVEL/3« 0 
manyhshanyh-i 
60  TO  1210 

1200  continue 

DIN()EX=-10  •  *AL0G10  ( TOJ/ (4t  *PI )  i 
IFUREV.NE.O)  RLEVEL=-10,»AuO610(TRL/(4.*PI) ) 

IFUNOEX  ,NE.  0)  WRITE  14*20)  DXNDEX 
IFUREV  .Nc.  0)  WRITE  (4r8D  R LEVEL 
60  TO  1Q0 

99«9  CALL  EXITG(UaTA) 

STOP 

ENO 
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XYZ 

SUBROUTINE  XYZ<IXTZ»XH»YH,2H,L»N»N»XSPAce,YSPACE»ZSPACE,LMN,X0,Y0, 
*20) 

C  *•**  A  SUBROUTINE  TO  OBTAIN  ARRAY  ELEMENTS  LOCATIONS 
DIMENSION  XM(1)»YHU)»2HU) 

200  format n 

LMNsL*M*N 

IF(lXY2  ,Nfe.  0)  00  TO  * 

IFUSPACf  ,LE«  0,0)  00  TO  3 
C  *9**  EQUALLY  SPACED  IN  X-OXRECTION 
1X0=0 

XO=FLOAT (L-i >  *XSPACE/2,0 

DO  1  1=1, L 

DO  11  JsliM 

DO  2l  K=1,N 

IXO=IXO+l 

XH ( I XO ) =FLOAT  U-l> *XSPACE«X0 
21  continue 
U  continue 
1  continue 

3  5F(YSPACE  ,LE.  0,0)  00  TO  5 

C  **••  tCl'ALLY  'fACEO  IN  Y-OlKECTlCN 
lYOsO 

YO=FLOAT (M-l ) *YSPACE/2,0 

DO  4  U=1»M 

00  14  Isl»L 

00  2*  K=1,n 

l.YO=IYO*l 

YH ( I YO) =FLOAT (M-l } *YSPACE-YO 
2*»  CONTINUE 
14  CONTINUE 

4  CONTINUE 

5  IFUSPACE  *LE.  0.0)  00  TO  6 

C  ****  EQUALLY  SPACED  IN  2-DIRECTION 
IZO=0 

20=FL0AT(N-1)*2S)  ACE/2,0 
00  7  K=1»N 
00  17  1=1 iL 
DO  2?  0=1, M 
120=120*1 

2H( 120) =FLOAT (K-l)*2SPACE*20 
27  CONTINUE 
17  CONTINUE 

7  CONTINUE 

6  return 

C  **•*  ELEMENTS  UNEQUALLY  SPACED 

8  CONTINUE 

xo=o.o 

YO=0»0 

20=0,0 

read  100,  UMN 
DO  9  I=1»LMN 

READ  100,  AN(I),YH(I) ,ZH(D 
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X0=X3*XHU> 

YOsyO+YHti) 

2^20+ZHU) 

9  CONTINUE 

XOSXO/FLOATCtHN) 
YOSYO/FlOAT(wHN) 
ZOSZO/FLOATaMN) 
DO  10  1=1 l CNN 

XH(I)=XHU)-XO 

YH(l)=tHU>-tO 

ZMU>SZH\l)-*0 

continue 

RETURN 

END 
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SUM 

SUBROUTINE  SUM(THETA»PHI»R2ERO»XH,Ytt,ZH,WLAHOA»LNN,OELA»R*RMAX#XMA 
*X#IsTEEK»SHAO) 

c  PERFORM  TRIPLE  summation 

DIMENSION  XH(l}»YHU)»ZH(i)'OELACl)tRU)»SHADU) 

DATA  RAO»PI/.17*5329E-01 .6.283185307/ 

TH=TH£TA*RAO 

PH=PHl*RAy 

ALPrtA=COS(PH)*SIN(TH) 

BETa=COS<PH) *COS (TM) 

OANmAs-SIN(PH) 

R£Al=0,0 

XlMsO.O 

IFUSTEER  «EQ.  0)  60  TO  ISA 
IF(RZERO  »6T.  0.0)  60  TO  60 
C  ****  PLANE  have  FRONT 
00  103  I=1.LHN 

AR6=PI* (XH( I ) *ALPHA«YH(1 )«B£TA«ZH (1 ) *6AMMA)/WLAmDA 
AR6sAR6*0ELA(l) 

REAu=REAL^bHAP ( I ) «R ( X ) *C0S ( AR6 ) /RMAX 
XlMsX 1M«SHA0 ( I ) «R ( X ) *SXN ( AR6 ) /RMAX 
103  CONTINUE 
60  TO  *1 

C  COMPUTE  STEERIN6  DIRECTIONS 

10*  CONTINUE 

00  105  1=1. LMN 

AK6uPl  *(  XH  U  )•  ALPHAS  YH  (X  )  «eETA*ZH  IX  )  *6AMMA  )/ WLAMDA 
OELA(X)sOELA(I)-AR6 
105  CONTINUE 

XMAx=FL0AT(LMN)**2 

ISTEERSI 

RETURN 

C  m*  SPHERICAL  *AV£  FRONT 

*0  continue 

alpha=aupha«rzero 
beTa=beta*»uero 
6ANmAS6AMMA*RZER0 
DO  203  1=1*  UMN 

AK6=  ( XH  ( I )  «ALPHA )  **24  (  YH  (  X  ) -BETA )  **£♦  (  ZH  (I )  "OAMmA  )  **2 

AR6=PI* (RZtKO-SQRT (AR6) ) /VLAMDA 

AR6sARG*t€LA(l) 

REA(.=REAL*SHAD(I)*K(I)*C0Si.Vl6)/RMAX 
XIMsXIM+SHAO t I ) *R ( l ) *SlNt AR6) / JUiAX 
203  CONTINUE 

*1  XMAx=REAL«KcAL+XIM*XXM 
RETURN 
END 
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c  beam  plot  -  rectangular 
c 

SUBROUTINE  BPLOTRUQEAN»ANGl£#!X»»NPTS,AFIX»AMIN,AmaX»ALONS»OLLUNG, 
•HOR.OATA) 

DIMENSION  ANGLE(l) 'DL(l) tHDft(l) 'DATA(i) 
data  RPl/310.2/ 
c 

C,  {BEAM  not  o:  vertical  pattern 

v  =  os  HORIZONTAL  PATTERN 

C  AFIXS  FIXED  ANGLE  (DEG) 

C  AMINS  MINIMUM  ANGLE  (DESS 
C  AMAXS  MAXIMUM  ANGLE  (DEG) 

C  ALONG!  LENGTH  OF  ANGLE  AXIS  (IN) 

C  OLLONG;  LENGTH  Of  LOSS  AXIS  (IN) 

C  HDftS  TITLE 
C 

XMIN  =  (G0»S.  -  AL0NG*RPI)/2. 

XMAx  =  XMIN  ♦  ALONG* RPI 

TMAx  =  3071.  -  ,5*RPI 

YMIN  =  YMAA  -  DLLONG*RPI 

CALL  0BUCTG(0ATA»XM1N*Y»IIN»XMAX»YMAX) 

OLMaX  s  io.*ullong 

call  SUBJEG(0ATA» AMlNfOLMAXf AMAX»0. ) 

ADEl  =  (AMAX-AMIN) /ALONG 
CALL  GRlOG(DATA»ASELf&0.«0»0) 

CALL  LAgELG(DATAr0'A0£U0»6.1) 

CALL  LAuELG(UATA»1«-I0«>0(3) 

CALl  TITLEG(0ATA»1x»*ANGLE  (DEG) *  »lb» *DB  DOAN  FROM  MRA**72»h0R) 

DO  20  I=lfNPTS 

IF (ANGLE ( l ) .LT • AMIN;  ANGLE(Z)  s  AMIN 
IF(ANGLEU).GT.AMAX)  INGLE (I)  a  AMaX 
IF(OL(I).LT.O«>  DL(i)  =  0. 

IF(i)L(I)  .GT.OLMAX)  BLU)  =  DLMAX 
20  CONTINUE 

CALL  SETSMGtOATA'SO'2,) 

CALL  LINESG(0ATA»NPTS,AN6LE»DL) 

CALl  SETSMtiuATA»30'l.) 

CALl  SETSMw(OATA»lG»I() 

Y  s  YMIN  •  l,25*KPI 
IF(IBEAm.EG.O)  GO  TO  SO 
X  s  (G09S,  -  21.*31.)/2. 

CALL  LEGNDG(0ATA»L*Y»21» ’VERTICAL  BEAM  PATTERN*) 

Y  =  T  -  .3*RPI 

X  s  (409b,  -  lB.*31.)/2. 

CALL  LEGNDG(OATA»X»Y»ia«*THETA  =  DEG*) 

X  s  X  ♦  31. *8. 

GO  TO  100 

90  X  s  (409b,  -  23.*3l.)/2. 

CALL  LEGN0G(0ATA»X»Yt23r*H0RIZ0NTAL  BEAM  PATTERN* ) 

Y  s  Y  -  *3*RPI 

X  s  (409b,  -  16.«31.)/2. 

call  LEGNOG(DATA»X»Y»10,«pHI  s  deg*) 

X  s  X  ♦  6,*3i» 

100  CALL  NUNBRG(0ATA»X»Y>0,1»AF1X) 

call  sctsmg(data»ig»o«; 
call  PAaEO(OATA»0,l#l) 

RETURN 

END 
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c  beam  plot  -  polar 
c 

SUBROUTINE  toPL0TP(IBEAM,AN6LE.DL*NPTS,AFIX*DLMAX»DLLCNG,H0R,DATA> 
DIMENSION  AN6LE (1) ,UL(1) *H0R(1) »0ATA(1) 

DATA  RP 1/510.7/ 

DATA  CEG/57 .2^5779/ 

I BEAM  NOT  05  VERTICAL  PATTERN 
=  01  HORIZONTAL  PATTERN 
AF1X;  FIXED  ANGLE  (DEG) 

DLMAx:  RANGE  Or  Dl  (OB) 

OLLONG:  RADIUS  (IN) 

hdr:  title 

XMlN  =  2047.5  -  OLLONG*RP| 

XMAX  5  XMIN  ♦  2.*DLL0NG*RFI 
YMAx  =  3071,  -  ,5*RPI 
YMIn  =  YMAX  -  2.*0LL0NG*RPI 
call  OBjCTG(OATA,XM1N,YMXN»XMAX,YMAX) 

CALL  SUB JEG ( 0 ATA . -OU.ON6 , -OlLONG  *  OLLONG , DLL0N6 ) 

CALL  RADIAL (OLLONG+ •  12  * DATA ) 

call  subueg (data, -olmax. -dlmax »dlmax»dlmax) 
call  SETSMG(OATA»30»2,) 
call  SETSMG(0ATA,63»1,) 
call  C2RARu(DATA,0. *0, ,DLMAx,0. ,360. ) 

CALL  SETSMG(0ATA,30fl,) 

00  10  1=0,170*10 
THETA  =  FL0AT(I)/D£G 
X  *  DLMAX*C0S( THETA) 

Y  S  DLMAX*S1N( THETA) 

10  CALL  SEGMTG(OATA»l»X»YfXt-Y) 

DO  20  1=1,10 
X  5  10. *1 

IF(X.GE.DLMAX)  GO  TO  30 
20  CALL  CIRAR»(OATA*0.,0.,OL«IAx-X*0,*360.) 

30  CALL  SETSMG(DATA,30*2.) 

THETA  =  (90.  -  ANGLE(1))/0EG 
R  =  DLMAX  -  OLU) 

X  =  R*COS(ThETA) 

Y  =  R»SIN(THETA) 

CALL  L1NESG(OATA*0*X*Y) 

00  40  I=2»nPTS 
THEtA  =  (90,-ANGLE(i))/OE« 

R  =  DLMAX  -  OL(!) 

X  =  R*COS( THETA) 

Y  =  R*SIN(IHETA) 

40  CALL  LINESG(DATA*1*X,Y) 

call  SETSMG(DATA*30»1,) 

CALu  SETSMG (DATA* 14*1,) 

Y  =  YMIN  -  160. 

X  =  (4095,  -  72.*3i.)/2. 

CALL  LEGN0G(DATA*X»Y*72*H0R) 
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Y  =  YMIh  •  1,25*RPI 
XF(IBEAH.Evi.O)  60  TO  50 
X  =  (4095,  -  21.-$31.)/2. 

CAM.  LE6ND6(DATA»X»Y»21»»VEkTXCAL  BEAM  PATTERN*) 

Y  s  Y  •  »3*RPX 

X  s  (4095,  -  18,*3l.>/2. 

CALL  LE6U06(DATA»X»Y.iB,»TH£TA  s  DEG') 

X  s  X  ♦  31. *8. 

60  TO  100 

SO  X  =  (4095,  -  23,«31,)/2. 

CALL  LC6N06(DATA'X*Y'25'*H0ftlZ0NTAL  BEAM  PATTERN* ) 

Y  s  Y  -  »3*KPX 

X  =  (4095,  -  X6,*31, )/2, 

CALL  LE6NCS(DATA#X»Y#16*  *RHl  s  0£6») 

X  s  X  ♦ 

100  CALL  NUMBR6(0ATA#X#Y»G*1»AFIXX 
XsRPI*DLL0N6/0LMAX*10 , 

YSYMAX425, 

00  120  XsiilO 
XX=lO*X 

IF(FLOAT(IX)  ,6E,  DLMAX)  SO  TO  ISO 
YsY-X 

120  CALL  NUMBR6(DATA»2032,5»Y»3,XX) 

130  CONTINUE 

CALL  SCTSMta(DATA;14»0,) 

CALL  PAGE6(0ATA»C,1#1) 

return 

end 
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SUBROUTINE  KAOIAU (RAD, DATA) 
DIMENSION  uATA(1) 

EXTERNAL  F0NT2 
DATA  0EG/S7.29578/ 

CALL  VEClS(DATArF0NT2f0) 

CALL  $ETSMG(OATA»Sl»i() 

CALl  NUMBR«(CATAf0.'RAO»l«0) 
0ANG2  =  l»/2»*31,/(RA0*27f*18) 
0ANC3  =  2»*DANS2 
00  1  KslOf 90f 10 
FPN  =  FLOAT (K) 

OANe  =  FPN/0E6  -  QANG2 
XO  s  RAD*SIN(DANS) 

TO  s  RAD*COS(DANG) 

THETA  3  -FPN 

CALL  SETSM0(0ATA»«6»THETA) 

CALL  NUMBRG(DATA'X0'Y0»2»K) 

1  CONTINUE 

DO  2  KsiOO* 260*10 
FPN  =  FLOAT (K) 

DANG  =  FPN/&E6  ♦  0AN63 
XO  =  RAD*S1N(0ANG) 

TO  =  RA0*C0S(DAN6) 

THETA  =  180.  -  FPN 
CALL  SETSM6(0ATAt46»THETA) 

CALL  NUmBR6(0ATA#X0#YQ#3»K) 

2  CONTINUE 

DO  3  K=270»350»10 

fpn  =  float (K) 

dang  =  FPN/DE6  -  DANGS 

XO  s  RA0«S1N(DANG) 

YO  =  RAD*CUS(DANG) 

THETA  s  360.  -  FPN 
CALL  SETSMG(0ATA'I»6'THETA) 

CALl  NUMBRu(0ATA»X0»Y0»3»K) 

3  CONTINUE 

CALL  SETSMG(DATA'51»0() 
call  SETSMG(DATA»«6»Ot) 
return 
end 


42 


SHADE 


SUBROUTINE  SHAD£(LMN'SKA0|XH'YH'2H'X0'Y0<Z0> 
DIMENSION  SHADU),XH(1)»YH(1)»ZH(1) 

PEDsO.S 

C  *•**  COSINE  SQUARE  SHADING  TECHNIQUE 
C AY 1SACOS ( SORT ( 0 a ) ) /XO 
c AY3=AC0S ( SURT ( o  a ) > /zo 
DO  99  I=1»LMN 

XA=  U . -PEO ) *COS (CAY 1*XH (I)) »*2+PE0 
XB=  (I  • -PCD)  *COS  (CAY9*ZH  ( I } )  «*24T>ED 
SHA0(X)=SHAU(I)*XA*XB 
99  CONTINUE 
RETURN 
ENO 


DELAY 

SUBROUTINE  DELAY (LMNrDELA) 

C  ***•  SELECTABLE  PHASE  DELAY  IN  RADIANS 
DIMENSION  UELA(i) 

100  FORMATO 

read  loot  (OELA(I)  asi»LMN) 

RETURN 

END 


RSPNSE 


SUBROUTINE  KS9NS£(LMNfR) 

C  RESERVED  TO  OBTAIN  DIFFERENT  RESPONSE  OF  DIFFERENT  ELEMENTS 
DIMENSION  HU) 

RETURN 

END 


REVB 


SUBROUTINE  KEVBiXt THiPH) 
Xs0.5 

return 

END 
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REVERSE  BLANK 


